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ABSTRACT 

o 
o 

Silicon nanoparticles (SNPs) have been proposed as the source of the observed "extended red 
emission" (ERE) from interstellar dust. We calculate the thermal emission expected from such 
particles, both in a reflection nebula such as NGC 2023 and in the diffuse interstellar medium 
(ISM). Pure neutral Si SNPs would emit at 16.4/im, while Si/Si0 2 SNPs (both neutral and 
charged) produce a feature at 20/im. Observational upper limits on the 16.4/^m and 20/im features 
in NGC 2023 impose upper limits of < 1.5ppm in pure Si SNPs, and or < 0.2ppm in Si/SiC>2 
SNPs. The observed ERE intensity from NGC 2023 then gives a lower bound on the required 
photoluminescence efficiency ??pl- For foreground extinction Ao.es^m = 1-2, we find i]pl > 5% for 
Si SNPs, or n PL > 24% for Si/Si0 2 SNPs in NGC 2023. Measurement of the R band extinction 
toward the ERE-emitting region could strengthen these lower limits. The ERE emissivity of the 
diffuse interstellar medium appears to require >42% (>33%) of solar Si abundance in Si/Si02 
(Si) SNPs. We predict IR emission spectra and show that DIRBE photometry appears to rule 
out such high abundances of free-flying SNPs in the diffuse ISM. If the ERE is due to SNPs, 
they must be either in clusters or attached to larger grains. Future observations by SIRTF will 
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O ' be even more sensitive to the presence of frcc-flying SNPs. 

Subject headings: dust, extinction — infrared: ISM: lines and bands — reflection nebulae: NGC 
2023 
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1. Introduction 

First detected in the Red Rectangle (Schmidt, Cohen, & Margon 1980), "extended red emission" (ERE) 
from interstellar dust consists of a broad, featureless emission band between ~5400 A and 9000 A, peaking at 
6100 < A p < 8200 A, and with a width 600 A < FWHM < 1000 A. The ERE has been seen in a wide variety 
of dusty environments: the diffuse interstellar medium (ISM) of our Galaxy, reflection nebulae, planetary 
nebulae, HII regions, and other galaxies (see Witt, Gordon, & Furton 1998 for a summary). The ERE 
is generally attributed to photoluminescence (PL) by some component of interstellar grains, powered by 
ultraviolet (UV) /visible photons. The photon conversion efficiency of the diffuse ISM has been determined 
to be near (10 ± 3)% (Gordon et al. 1998; Szomoru & Guhathakurta 1998) assuming that all UV/visible 
photons absorbed by interstellar grains are absorbed by the ERE carrier. The actual photoluminescence 
efficiency ?7pl of the ERE carrier must exceed ~ 10%, since the ERE carrier cannot be the only UV/visible 
photon absorber. 

Various forms of carbonaceous materials - hydrogenated amorphous carbon (HAC) (Duley 1985; Witt & 
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Schild 1988), polycyclic aromatic hydrocarbons (PAHs) (d'Hendecourt et al. 1986) 1 quenched carbonaceous 
composite (QCC) (Sakata et al. 1992), C 60 (Webster 1993), coal (Papoular et al. 1996), PAH clusters 
(Allamandola, private communication), and carbon nanoparticles (Seahra & Duley 1999) 2 - have been 
proposed as carriers ol ERE. However, most candidates appear to be unable to simultaneously match the 
observed ERE spectra and the required PL efficiency (see Witt et al. 1998 for details). 

Very recently, Witt et al. (1998) and Ledoux et al. (1998) suggested crystalline silicon nanoparticles 
(SNPs) with 15A - 50A diameters as the carrier on the basis of experimental data showing that SNPs 
could provide a close match to the observed ERE spectra and satisfy the quantum efficiency requirement. 
It was estimated by Witt et al. (1998) and Ledoux et al. (1998) that SNPs account for < 5% of the total 
interstellar dust mass, with Si/H w 6 ppm. Zubko et al. (1999) modelled the interstellar extinction curve 
taking SNPs as an interstellar dust component containing Si/H w 18 ppm, or ^50% of the solar abundance 
[Si/H] Q = 36 ppm. More recently, Smith & Witt (2001) have further developed the SNP model for the ERE, 
concluding that the observed ERE in the diffuse ISM can be explained with Si/H = 6 ppm in SiCVcoated 
SNPs with Si core radii a w 17.5 A. 

The purpose of this paper is to test the SNP hypothesis. We calculate the IR spectra for SNPs in a 
reflection nebula - NGC 2023 - and compare with observations. We show that the SNPs in NGC 2023 
contain Si/H < 1.5 ppm. We re-estimate the minimum Si depiction in SNPs required to account for the 
observed ERE intensity in the diffuse ISM, and calculate the IR emission expected from such particles. We 
show that existing DIRBE photometry appears to rule out the abundances of free-flying SNPs required 
to account for the ERE emissivity of the diffuse ISM. Future observations by the Space Infrared Telescope 
Facility (SIRTF) will be even more sensitive to the abundance of SNPs in the diffuse ISM. 

In §2 we discuss the optical properties of silicon nano-crystals and glassy Si02 (§2.1) as well as the 
heat capacities of pure Si and Si core-SiC>2 mantle grains (§2.2). In §3 we carry out calculations for the IR 
emission spectra of SNPs (with or without SiC>2 coatings), and discuss their implications. In §4 we discuss 
the effects of grain shape, mantle thickness, and sources of uncertainties. Our conclusions are presented in 
§5- 

2. Grain Physics 

Experimental studies indicate that high photoluminescencc efficiencies are only observed from Si nanocrys- 
tals when their surfaces have been "passivated" by oxidation or hydrogenation - otherwise electron-hole pairs 
recombine nonradiatively (see, e.g., Kovalev et al. 1999). Accordingly, SNPs with oxide coatings are likely 



1 High photolumincsccnce efficiencies can be obtained by PAHs. Arguments against PAHs as ERE carriers include: (1) 
presence of sharp structures in the luminescence spectra of individual PAH molecules in contrast to the featureless nature of the 
interstellar ERE spectra; (2) lack of spatial correlation between the ERE and the PAH IR emission bands in the compact HII 
region Sh 152 (Darbon et al. 2000) and in the Orion Nebula (Perrin & Sivan 1992); (3) ERE detection in the Bubble Nebula 
where no PAH emission has been detected (Sivan & Perrin 1993); (4) nondctcction of ERE emission in reflection nebulae 
illuminated by stars with effective temperatures T* < 7000 K (Darbon, Perrin, & Sivan 1999) whereas PAHs emission bands 
have been seen in such regions (e.g., see Uchida, Sellgren, & Werner 1998) and are expected for the PAH emission model (A. 
Li & B.T. Drainc 2001, in preparation). Argument (1) may not be fatal since a featureless band may result from a mixture of 
many individual PAH molecules and ions. 

2 Seahra & Duley (1999) argued that small carbon clusters were able to meet both the ERE profile and the PL efficiency 
requirements. However, this hypothesis appears to be ruled out by non-detection in NGC 7023 of the 1/im ERE peak (Gordon 
et al. 2000) predicted by the carbon nanoparticle model. 
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to be of primary interest as the source of the ERE. However, we will also discuss the infrared emission of 
pure Si SNPs in order to show that even if the Si core dominates the IR emission, the IR spectral signature 
is still conspicuous. 



The optical properties of crystalline Si nanoparticles are controversial. Some studies (see Yoffe 2001 for 
a recent review) have concluded that the optical properties of Si nanoparticles differ substantially from those 
of bulk Si. This has been attributed to the quantum confinement effect (Wang & Zunger 1994; Tsu, Babic, 
& Ioriatti 1997). Koshida et al. (1993; hereafter KKS93) have published dielectric constants for crystalline 
nanosilicon which are considerably smaller than those of bulk-Si. We find that the absorption and reflectivity 
measurements for porous silicon (e.g., Koshida et al. 1993, de Filippo et al. 2000) appear to be inconsistent 
with the KKS93 optical constants; instead, they can be approximately reproduced using a mixture of Si and 
voids, if the Si component is described using the optical constants of bulk Si, as previously found by Kovalev 
et al. (1996), Theiss (1997), Lerondel, Madeore, & Muller (2000), and Diesinger, Bsiesy, & Herino (2001). 
Oxide layers are generally also present in laboratory samples of porous Si; we find that the absorption and 
reflectance spectra can also be satisfactorily reproduced using Bruggeman effective medium theory (Bohren 
& Huffman 1984) for a mixture of voids, Si0 2 , and bulk-Si. 

We take the following "synthetic" approach to obtain the complex refractive index to(A) = m' + im" 
for bulk Si at low temperatures, with no electron-hole pairs present (the contribution of electron-hole pairs 
will be treated separately). For 0.01 < A < 1.2/im we take m" of Adachi (1999) for bulk crystalline Si. For 
1.2 < A < 3.6/mi we set to" = 2.4 x 10~ 8 estimated from the (room temperature) absorption coefficient 
a w 10 _3 cm _1 of bulk Si between 1 and 3^m (Gray 1972). For 3.6 < A < 25/xm we take those of Palik 
(1985) for crystalline bulk Si; extrapolation is then made for A > 25/um. 3 After smoothly joining the adopted 
to", we calculate to' from to" through the Kramers-Kronig relation (Bohren & Huffman 1983, p. 28). 

The Si long-wavelength absorption properties will depend on whether any free electrons or holes are 
present. SNPs arc thought to luminesce only when uncharged and containing exactly one electron-hole 
pair, which recombines radiatively in ~ 10~ 3 s (Smith & Witt 2001). The thermal energy remaining after 
luminescence will then be radiated by a SNP with no free electrons or holes, at least in the diffuse ISM where 
the cooling time is short compared to the time between absorption of UV photons. For these SNPs, we use 
the optical properties of bulk Si. 

Charged SNPs will contain at least one free electron or hole, which we assume contributes to the 
dielectric function 



3 We do not adopt those of Palik (1985) for A > 25/im since they were measured at T k» 300 K, and the A > absorption is 

dominated by thermally-excited electron-hole pairs. In the absence of such pairs, we approximate m"(A) « m,"(25^tm)(25/im/A). 
The IR bands of Si nanoparticles may become stronger due to the symmetry-breaking surfaces (e.g., see Hofmeister, Rosen, 
& Speck 2000). New IR bands, forbidden in bulk Si, may also appear (Adachi, private communication). Unfortunately, IR 
measurements of Si nanoparticles are unavailable, so we must use the optical properties of bulk Si. 



2.1. Optical Properties 
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where we adopt an effective mass m e ff ~ 0.2m e for electrons or holes in Si. 4 Scattering by the Si boundary 
will result in r « cl/vf, where we take vf ~ 10 s cms" 1 as the typical velocity 5 

In Figure 1 we plot the resulting optical constants for crystalline SNPs, both for neutral SNPs and 
SNPs with a radius a = 10 A and Z = +1. We see that the hole dominates the absorption for A > 1/xm, 
and completely overwhelms the weak vibrational bands in the 7-25/xm region which provide the infrared 
absorption in neutral Si. For comparison, we also show the optical constants from KKS93. 

For glassy SiC>2, we take m" from Palik (1985) for 0.01 < A < 7/im but with one modification: for 
0.15 < A < 3.6/Ltm we set m" = 1.0 x 10 -4 based on the absorption coefficient measured by Harrington et al. 
(1978). For 7 < A < 500/xm we take m" from Henning & Mutschke (1997) for glassy Si0 2 sample at 100 K 
(the optical properties of glassy SiC>2 are not sensitive to temperature, see Henning & Mutschke 1997). For 
A > 500^m, we approximate m"(A) « m"(500^m)(500/im/A). Again, the real part m' is calculated from the 
Kramers-Kronig relation. The results are also presented in Figure 1 . 

In Figure 2 we show the absorption cross sections per unit volume calculated for both spherical and 
spheroidal (a) pure Si grains, and (b) Si/Si02 grains, with the Si in a confocal spheroidal core containing a 
fraction (0.8) 3 = 51% of the volume. We also plot in Figure 2c the radiation fields for the diffuse ISM and 
for NGC 2023 (60" S of HD 37903; see §3.1). 

Figures 2a and 2b clearly show that in the wavelength range where Si or Si/Si02 grains absorb most (A < 
0.25/xm) the grain shape can substantially affect the wavelength-dependence of the ultraviolet absorption. 
However, we will see below (§4) that the integrated energy absorption rate, and the resulting emission 
spectrum, is surprisingly insensitive to grain shape. 

2.2. Enthalpies 

The experimental specific heat of bulk crystalline Si can be approximated by a Debye model (with 
dimensionality n = 3) and Debye temperature 9 = 530 K, while SiC>2 glass can be approximated by a model 
where 1/3 of the vibrational modes are distributed according to a Debye model with 6 = 275 K, and 2/3 of 
the modes according to a Debye model with 6 = 1200 K. 

For very small particles at low temperatures, the discrete nature of the vibrational spectrum becomes 
important. Let N^ om and be the number of atoms in the Si core and Si0 2 mantle, and JV atom = 

^atom + ^atom be the total number of atoms. For pure Si clusters, we set N m = 3N^ om — 6 vibrational modes 
where the "—6" term allows for the translational and rotational degrees of freedom since the energy from 
photon absorption is only distributed among the vibrational modes. For Si/Si0 2 nanoparticles, we assume 
N m = SN^ om — 3 modes distributed according to a 9 = 530 K Debye model, N m = 2iV at 1 Jj — 2 modes 
distributed according to a 9 = 1200 K Debye model, and N m = N a t £ — 1 modes distributed according to a 
9 = 275 K Debye model. We assume the mode frequencies to be distributed following eqs.(4-6) and eq.(ll) 
of Draine & Li (2001). The specific heat is calculated treating the modes as harmonic oscillators, with the 



4 For a positively charged SNP, the relevant mass is the effective mass of a hole. There are two degenerate valence band 
maxima, with hole masses 0.49m e and 0.16m e (Ashcroft & Mermin 1976, p. 569); we adopt 0.2m e as a representative value. 
The conduction band minima are characterized by electron effective masses of 1.0 and 0.2m e ; these would be relevant for 
negatively-charged SNPs. 

5 Our estimate of r is by no means exact, but the value of r is unimportant in the case of oxide-coated SNPs, since the oxide 
layer dominates the IR emission. As mentioned above, interstellar SNPs are thought to be oxide-coated. 
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Fig. 1. — Optical constants m! (upper panel), m" (lower panel) of neutral Si (solid lines), charged Si (Z = +1, a = 10 A; 
dashed lines), and SiC>2 glass (dot-dashed lines). Also plotted are those of Koshida et al. 1993 (dotted lines; labelled by KKS93). 

continuum limit used for large particles and high energies (see Draine & Li 2001). 

3. IR Emission Spectrum and Its Implications 

Let Zsi be the amount of Si in the grains (by number) relative to total hydrogen. The IR emissivity per 
unit solid angle per H nucleon from the SNPs is 

^ = TT74 C abs(a,A) / dT B X (T) P(a,T) (3) 
Nsi{a) Jo 

where iVg;(a) is the number of Si atoms in a grain of radius a [7Vgj = (7ra 3 /21mH)psi f° r pure Si nanoparticles, 
and -ZVgi = Tra^psi/^lma + 7r(a 3 — agi)psi0 2 /45mH for grains with Si cores of radius asi and SiC>2 mantles, 
where psi, Psi0 2 are the mass densities of crystalline Si (« 2.42gcm -3 ) and glassy SiC>2 (~ 2.2gcm~ 3 ), 
respectively, and mu is the mass of a hydrogen atom]; C a bs(a, A) is the absorption cross section of the 
spheroidal core-mantle grain at wavelength A, calculated using the optical constants discussed in §2.1; B\(T) 
is the Planck function; P(a, T)dT is the probability that the grain temperature will be in [T, T + dT] . For a 
given radiation field, we calculate P(a, T) for small grains employing the "thermal-discrete" method (Draine 
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Fig. 2. — (a): Absorption cross sections per unit volume of spherical Si grains (solid line), 5:1 prolate Si grains (dashed line), 
5:1 oblate Si grains (dotted line). Also plotted are the spherical, 5:1 prolate, 5:1 oblate grains using the optical constants of 
Koshida et al. (1993) (labelled by KKS93). (b) Absorption cross sections per unit volume of spherical Si/Si02 grains (solid 
line), 5:1 prolate Si/Si02 grains (dashed line), 5:1 oblate Si/SiC>2 grains (dotted line). All grains are taken to have the same Si 
core volume fraction 51% (for the spherical case, agj = 0.8a). (c): The solar neighbourhood interstellar radiation field (Mathis, 
Mezger, & Panagia 1983) and the NGC 2023 radiation field at a position 60" S of HD 37903 (see §3.1). 



& Li 2001). 6 We characterize the intensity of the illuminating starlight by Xi the intensity at 1000 A relative 
to the Habing (1968) radiation field. For NGC 2023 the spectrum is assumed to be a 22000 K dilute blackbody 
cutoff at the Lyman edge; for the diffuse interstellar medium we take the spectrum of Mathis, Mezger, & 
Panagia (1983), with \ = 1-23. 



3.1. NGC 2023 

The reflection nebula NGC 2023, at a distance D = 450pc, is illuminated by the B1.5V star HD 37903 
(T e g = 22,000K, L+ = 7600i Q ). A modest HII region surrounds the star, beyond which is a photodissoci- 



6 To be precise, the energy released in the form of PL is not available as heat. In NGC 2023, with a mean energy (/ivJere ~ 
1.8eV for ERE photons, the fraction of absorbed photon energy lost to photoluminescence i)pL('« / )ERE/('"')abs i s only 0.21»;pl 
for pure Si ((/if) a bs = 8.6 eV) or 0.19t;pl for Si02-coated Si ((/i^) a bs = 9.3 eV). Even for ijpl — ► 1, this is a minor correction 
which does not alter our conclusions. 
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ation region at an estimated distance ~ 5 x 10 17 cm from the star. The radiation field is expected to have 
an intensity x ~ 5000 at this distance from the star. While x ~ 5000 is consistent with models to reproduce 
the observed IR and far-red emission from UV-pumpcd H2 in the bright "emission bar" 80" S of HD 37903 
(Draine & Bertoldi 1996, 2000), the integrated infrared 5-60/xm surface brightness at a position 60" S of 
HD 37903 (Figure 5) points to a lower value of \- We will assume that the dust producing the measured 
infrared spectrum and ERE is illuminated by a radiation field with x ~ 2000. 7 

In Figures 3 and 4 we show temperature distribution functions P(a, T) for pure Si nanoparticlcs and 
Si/SiC>2 nanoparticles of various sizes. In each case we show P(a, T) for both neutral and charged grains. 
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Fig. 3. — Temperature distribution functions for pure Si spheres of various radii, using the optical properties of (a) neutral 
particles and (b) particles with charge Z = +1, exposed to the radiation field in NGC 2023. Curves are labelled by grain radius 
a. 

The 2.4-45/im spectrum has been obtained by Verstraete et al. (2001) using the Infrared Space Observa- 
tory, using apertures ranging from 14" x 20" (for 2.4-12/zm) to 20" x 33" (for 29-45//m), at a position 60" S 
of HD 37903. Far-infrared photometry with a 37" beam centered 60" S of HD 37903 has been obtained from 
the Kuiper Airborne Observatory by Harvey et al. (1980). The observed IR spectrum is shown in Figure 5. 
Proposed ERE carriers must not produce IR emission in excess of what is observed. 



The PDR in NGC 2023 is optically thick to the illuminating starlight. Allowing for forward scattering, 



7 As explained in §4, our analysis is insensitive to the precise value assumed for \. 
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Fig. 4. — Temperature distribution functions for Si/Si02 spheres of various radii, with the Si cores assumed to be (a) neutral 
and (b) charged with Z = +1, for grains exposed to the radiation field in NGC 2023. Curves are labelled by the Si core radius 
a Si (~ 0.8a where a is the Si/Si02 radius). 



we estimate that the effective attenuation cross section per H nucleon for 3-11 eV photons is ~ 0.5 x 
10 -21 cm 2 H _1 , and we approximate the PDR as a slab of optically-thin dust and gas with Nn rs 2x 10 21 cm~ 2 . 
We further assume a limb- brightening factor 1/ cos (9 rj 2, where 6 is the angle between the slab normal and 
our line-of-sight (our final results are insensitive to the adopted values of Nn and 1/ cost?; see §4). With the 
above assumptions, j\, the power radiated per H nucleon per unit solid angle per unit wavelength A (see 
eq. [3] ) is related to the observed intensity I\ by 

h = H* jxlQ -0.*Ai „ 4 x 1Q 21 cm -2 jxW -0AA, (4) 

cos a 

where A\ is the extinction between the point of emission and the observer. 

The foreground extinction A\ is uncertain. The line-of-sight to HD 37903 has E(B — V) w 0.35 and 
R v = A V /E(B - V) = 4.1 (Cardelli, Clayton, & Mathis 1989), corresponding to A . 68fJ , m ~ 1.2. For the 
emission bar 80" S of HD 37903, an extinction j4o.68/xm ~ 3.2 has been determined from the relative strengths 
of K band and far- red H2 emission lines (Draine & Bertoldi 2000). For purposes of discussion, we will assume 
that the extinction to the reflection region 60" S of HD 37903 (where the ERE has been measured) is likely 
to be in the range 1.2 < Ao.^m 3.2. 

To estimate 7 ERE (= J lf RE d\) 60" S of HD37903, we take (1) the measured surface brightness profile 
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from Witt, Schild, & Kraiman (1984, Figure 7); (2) h.m^l h.^m = W ' 06 Kes^J K(>2nm = 0.91 80" 
S of HD37903, where F* is the observed stellar flux (Witt & Schild 1988, Table 1); (3) the ERE spectral 
profiles 62" and 84" ENE of HD37903 (Witt & Boroson 1990, Figure 4). 8 From this we estimate 7 ERE w 
9 x 10~ 5 erg cm" 2 s" 1 sr" 1 60" S of HD37903. The ERE emissivity per H nucleoli is 

p f fER.EinO.4ylA J\ 

•ere / .ERE dA = J a 1U « 2.3 x l(T 26 10 - 4A °— erg s" 1 sr-^" 1 . (5) 

Proposed ERE carriers must provide this level of emission. 




A (/xm) 

Fig. 5. — Predicted IR emission spectrum for Si/H=12 ppm in (a) neutral and (b) Z = +1 pure Si nanoparticles in NGC 2023. 
For the neutral Si nanoparticles, the emission is nearly independent of particle size; for the charged Si nanoparticles, spectra 
are labelled by grain radius. Heavy solid line: ISO SWS spectrum (Verstraete et al. 2001). Stars: KAO photometry (Harvey 
et al. 1980). For the 16.4/^m feature not to exceed the ISO spectrum, the upper limit on neutral Si SNPs is Si/H < 0.8 ppm. 
The upper limit on charged Si SNPs alone is Si/H < 1.5 ppm. The upper limit on Si SNPs (both neutral and charged) is Si/H 
< 1.5 ppm. 



The ERE spectrum for 60" S peaks at 6800 A (Witt & Schild 1988) as do the spectra 62" and 84" ENE of HD37903. 
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3.1.1. Si Nanoparticles 

We first calculate the IR emission spectra for crystalline Si grains illuminated by radiation from HD 
37903 with intensity \ — 2000. As shown in Figure 3, the radiation field is so strong that neutral Si grains 
with a > 30 A attain a steady-state temperature « 820 K, while a > 30 A charged Si grains are heated to 
< 200 K. The single-photon heating effect (Greenberg 1968) is evident for a^i — 10 A. In Figure 5 we present 
the predicted spectra for Zsi = 12 ppm in Si nanoparticles 9 with radii a = 10, 20, 30, 40, 50 A. Grains with 
a < 100A are in the Raylcigh limit for A > 912A, and thus their "equilibrium" temperatures are independent 
of a; therefore, for a given nanoparticle abundance Zsi, the resulting IR intensity is almost independent of 
grain size. 

Although the fundamental lattice vibration of crystalline Si has no dipole moment and thus is IR 
inactive (i.e., there is no interaction between a single phonon and IR radiation), multi-phonon processes 
produce prominent bands at 6.91, 7.03, 7.68, 8.9, 11.2, 13.5, 14.5, 16.4, and 17.9^m (Johnson 1959). As 
seen in Figure 5a, for Zsi = 12 ppm the 16.4^m emission feature is calculated to be about 15 times stronger 
than the observed spectrum if the grains remain neutral. In order not to exceed the observed spectrum near 
16.4^m, NGC 2023 must have Z si < 0.8 ppm in neutral Si SNPs - less than 1/15 of the abundance of SNPs 
to explain the ERE in the diffuse ISM. If, however, the pure Si particles have charge +1, then the thermal 
emission from the hole results in an emission continuum peaking near 15-20/zm. We see from Figure 5 that 
an abundance Zsi = 12 ppm in charged SNPs would result in emission which would exceed the observed 
spectrum by about a factor ~ 10 near 20^m. Thus we obtain an upper limit Zsi < 1.5 ppm in pure Si SNPs 
(all charge states) in NGC 2023. 

Let ?7pl be the photoluminescence efficiency, 7si be the UV/ visual photon absorption rate per Si (in 
the 912-5500A wavelength range), and (/i^)ere be the mean energy of ERE photons. Since the NGC 
2023 ERE peaks at - 6800A (Witt & Boroson 1990), we take (/m/) E re « 1.8cV. Illuminated by starlight 
of T cff =22,000K, we calculate 73; w 7.43 x 10~ 6 (x/2000) s" 1 Si" 1 for neutral Si nanoparticles. The ERE 
cmissivity per H is then j' ERE = (ft.z^)ERE7si?7PLZsi/47r w 1.70 x 10~ 18 t^plZs; erg s _1 sr _1 H _1 . From 
eq. (5) we thus require r?pLZsi ~ 1-35 x 10~ 8 10 a4A ° - 68 e m . Correspondingly, for charged Si particles, we 
obtain 7si w 7.04 x 10- 6 (x/2000) s" 1 Si" 1 , j ERE « 1.62 x 10~ 18 ry PL Z S i erg s" 1 sr" 1 H" 1 , and 7?p L Zsi ~ 
1.42 x io- 8 10°' 4Ao 68 *"". Results are summarized in Table 1. 

Since we have seen that Zsi < 0.8 ppm, it follows that pure Si neutral nanoparticles could produce the 
observed ERE only if the photoluminescence efficiency ?7pl > 0.017 x I0 - 4j4 ° 68 ^ m . For charged Si grains 
(Z S i < 1.5 ppm), r] PL > 0.010 x l0 0AA o w is required. However, according to Smith & Witt (2001), those 
charged grains will not emit ERE, i.e., 77PL ~ 0. 

3.1.2. S1/S1O2 Nanoparticles 

It seems unlikely that pure Si nanoparticles would exist in interstellar space without some degree of 
oxidation. Si nanocrystals in air spontaneously form an overlayer of Si02 which is usually <~ 10A thick 
(Ledoux et al. 1998). Witt et al. (1998) propose that SNPs in the form of Si core-SiC>2 mantle structure are 
a product of the initial dust formation process in oxygen-rich stellar mass outflows. Laboratory studies show 



9 As discussed in §3.2, we estimate that the observed ERE from the diffuse ISM would require Zgi > 12 ppm in Si SNPs, or 
Z Si > 15 ppm in Si/SiQ 2 SNPs. 
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Fig. 6. — As in Figure 5, but for Si/Si02 SNPs, for various values of the Si core radius asi = 0.8a. Upper panel (a) — Model 
spectra calculated for neutral Si/SiC>2 grains containing Si/H=15 ppm, with radii a$i = 10, 20, 30, 40, 50 A (see text). Spectra 
for asi = 20, 30, 40, 50A are nearly indistinguishable. Lower panel (b) - As in (a), but for charged Si/Si02 grains. The ripple 
features at A > 30/im are caused by the Si02 optical constants adopted from Henning & Mutschke (1998). 

that the oxidation of the Si core is a self-limiting process and, in the final stage of oxidation, the thickness 
of the oxide layer represents about 10% of the total particle diameter (Ledoux et al. 2000). 

We now calculate the IR emission spectra for Si core-Si02 mantle particles (both neutral and charged). 
Following Ledoux et al. (2000), we assume asi = 0.8a where asi and a are respectively the radii of the Si 
core and the Si/SiC>2 entire grain. In Figure 6 we show theoretical spectra for agi — 10, 20, 30, 40, 50A with 
Zsi = 15 ppm (estimated in §3.2 for Si/SiC>2 grains to account for the ERE in the diffuse ISM) together with 
the observed spectrum for NGC 2023. As illustrated in Figure 6 the Si-0 modes produce strong features 
at 20/im, 12.5/im and 9.1^im. The most conspicuous feature is the strong and broad 20/im band which, for 
Zsi = 15 ppm, is ^80 times stronger than the observational data. To depress the 20/xm emission feature 
to a level not in contradiction with the observational spectrum, NGC 2023 must have Zsi < 0.2 ppm in 
SiCVcoated Si grains (all charge states). 

The 912-5500A photon absorption rate for these Si/Si0 2 neutral SNPs is 7si « 6.19xl0~ 6 (x/2000) s^ 1 Si -1 . 
The ERE emissivity per H is then (/if }ERE7Si?7PLZsi/47r w 1.42 x 10~ 18 rypLZsi erg s -1 sr _1 H _1 . From 
cq. (5) we thus require ^p L Z Si w 1.62 x 10~ 8 10 a4A ° 68 ^ m . Since Z Si < 0.2 ppm the observed ERE in- 
tensity in NGC 2023 (eq. [5]) requires a photoluminescence efficiency rjph > 0.081 x l0 - 4j4 ° 68 f m if due to 
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Si corc-Si02 mantle neutral nanoparticles. Similarly, we obtain for Si/SiC>2 charged SNPs 7s; ~ 6.00 x 
10- 6 (x/2000) s- 1 Si" 1 , j ERE w 1.37 x 10~ 18 r7P L Z S i erg s" 1 sr" 1 H" 1 , r/ PL Z si w 1.68 x 10- 8 10°- 4A ° - BS ^, 
and ?7pl > 0.084 x 10°- 4A ° ■ 68 /"». Again, according to Smith & Witt (2001), those charged grains will have 
i] PL « 0. Results are also summarized in Table 1. 



3.2. Diffuse Interstellar Medium 

Using Pioneer 10 and 11 observations, Gordon et al. (1998) estimated the ERE cmissivity per H nucleon 
of the high Galactic latitude (HGL) diffuse ISM to be j ERE « 1.4 x 10~ 26 ergs" 1 sr" 1 !!" 1 . Szomoru & 
Guhathakurta (1998) obtained (J ERE ) = 1.2 x 10~ 5 ergs" 1 sr^cm" 2 toward 3 cirrus clouds with (N H ) = 
7.3 x 10 20 cm~ 2 , giving j ERE « 1.6 x 10 _26 ergs _1 sr _1 H _1 . Evidently the ERE carrier must provide an 
emissivity j ERE w 1.5 x 10 -26 erg s" 1 sr _1 H _1 in the diffuse ISM. The ERE appears to peak at -6000 A 
(Szomoru & Guhathakurta 1998); if due to SNPs, this would indicate 2a « 28 A (Ledoux et al. 2000), or 
a si w 0.8a « 11 A. 



3.2.1. Si Nanoparticles 

The 912-5500A photon absorption rates 7s; for pure Si neutral and charged grains are « 4.94 x 
10- 9 (x/l-23)s- 1 Si" 1 , and«4.68xl0- 9 (x/1.23)s~ 1 Si" 1 , respectively. Thus j ERE « 1.32x 10- 21 (x/1.23)r7 PL Z S i erg s" 1 sr" 1 
and j ERE « 1.25 x 10~ 21 (x/l-23)r/PLZsi erg s _1 sr _1 H _1 , respectively, for neutral and charged Si grains. 
Therefore, since t^pl < 100%, we require Zsi > 12 ppm (> 33% of [Si/H] Q ) for both neutral and charged 
Si grains. In Figure 7 we plot the temperature probability distribution functions for these grains, and IR 
emission spectra are shown in Figure 8. 

The DIRBE 25/xm photometry is inconsistent with Zsi ~ 12 ppm in pure Si neutral nanoparticles, with 
an upper limit Zsi < 4 ppm. For charged Si grains, the predicted IR spectra arc not in conflict with the 
DIRBE photometry. However, we note again that, according to Smith & Witt (2001), those charged grains 
will not emit ERE. Results are summarized in Table 1. 



3.2.2. Si/Si0 2 Nanoparticles 

Similarly, we have also considered Si core-Si02 mantle grains with a s ; = 0.8a, assuming the Si core to 
be either neutral or have charge +1. 

The 912-5500A photon absorption rates 7 S ; for Si/Si0 2 neutral and charged grains are w 3.87 x 
10- 9 (x/l-23)s- 1 Si" 1 , and« 3.77xl0- 9 (x/l-23)s- 1 Si" 1 , respectively. Thus j ERE « 1.04x 10- 21 (x/1.23)r/ PL Z si erg s" 1 sr~ ] 
and j ERE w 1.00 x 10~ 21 (x/l-23)?7PLZsi erg s _1 sr _1 H^ 1 , respectively, for neutral and charged Si/SiC>2 
grains. Therefore, since t?pl < 100%, we require Zsi > 15 ppm (> 42% of [Si/H] Q ) for Si/Si02 grains. In 
Figure 9 we plot the temperature probability distribution functions for these grains, and IR emission spectra 
are shown in Figure 10. Results are also summarized in Table 1. 

Evidently, the DIRBE photometry rules out Zsi ~ 15 ppm in Si/Si02 nanoparticles (both neutral and 
charged) with radii asi < 35 A. The DIRBE photometry places an upper limit Zsi < 2 ppm on the abundance 
of the asi ~ 11 A particle size which appears to be required by the observed wavelength of peak emission in 
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Fig. 7. — Temperature distribution functions for pure Si grains illuminated by the local interstellar radiation field (Mathis, 
Mezger, & Panagia 1983). Upper panel: neutral Si grains. Lowel panel: charged Si grains. For 10A (upper panel), 10— 50A 
(lower panel) a dot indicates the first vibrationally-excited state (see Draine & Li 2001). 

the diffuse ISM. Thus, we conclude that free-flying SNPs can account for no more than 13% of the observed 
ERE from the diffuse ISM. 

SIRTF will be able to obtain diffuse interstellar cloud emission spectra which should either detect the 
16.4 or 20^jm spectral features, or place much stronger constraints on the abundance of SNPs. 

4. Discussion 

Si/SiC>2 core-mantle grains are heated primarily by absorption by the Si core (the photoluminescence of 
which is responsible for the ERE), since the SiC>2 mantle is almost nonabsorptive for 0.12 < A < 4/mi; the 
core-mantle grains are cooled primarily by the SiC>2 mantle, since m" is so small for Si in the infrared (see 
Figure 1). Therefore, increasing the amount of Si02 present (while holding the pure Si cores constant) will 
lead to only small changes in the IR emission resulting from the small decrease in average grain temperature 
(e.g., Si/SiC>2 grains with equal numbers of Si atoms in the core and in the mantle [asi ~ 0.668a] in NGC 
2023 [x=2000] have a steady-state temperature ~ 85K in comparison with «96K for grains with asi = 0.8a). 
The limits on rjp^ will be only slightly affected by increasing the thickness of the SiC>2 mantles. 
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A (/xm) 

Fig. 8. — Predicted IR emission spectra for pure Si SNPs with Zsi = 12 ppm in the diffuse ISM, for (a) neutral Si grains, and 
(b) Z = +1 Si grains. Curves are labelled by grain radius. Also plotted are the DIRBE photometry (diamonds; Arendt et al. 
1998) and the HGL model spectrum of Li k, Draine (2001). 

In §3 the grains are modelled as spherical. To investigate the sensitivity to shape variations, we have 
also calculated the IR emission spectra for (1) 5:1 prolate grains; (2) 5:1 oblate grains; (3) grains with a 
distribution of spheroidal shapes with dP/dL\\ — 12L|| [1 — L||] 2 (Ossenkopf, Hcnning, & Mathis 1992) where 
L\\ is the so-called "depolarization factor" parallel to the grain symmetry axis (for spheres L|| = l/3); this 
shape distribution peaks at spheres and drops to zero for the extreme cases Ln — > (infinitely thin needles) or 
Lii — > 1 (infinitely flattened pancake). For core-mantle grains, we assume confocal geometry, with the above 
dP/dL\\ applying to the outer surface. For illustration, we present in Figure 11 the IR spectra calculated 
for neutral Si/SiC>2 grains illuminated by the NGC 2023 radiation field. It is clear that the shape effects are 
minor: the major change is a slight shift of the peak wavelength of the 20^m Si-0 feature. 

In a reflection nebula such as NGC 2023, the IR surface brightness due to SNPs J I\dX oc A^hxZs;/ cos 8. 
While the fraction of the IR power radiated in the 20/xm Si02 feature (or the 16.4/zm Si feature in the case 
of pure Si nanoparticles) depends on \ (through the grain temperature), this dependence is relatively weak 
for modest changes in x- Therefore the upper limits on the 16.4 and 20/im features in NGC 2023 give upper 
limits on A^nX^Si/ cos#; these limits, plus the observed ERE oc 77PL-^HXZsi/ cos^ then give lower limits on 
t]pl, independent of the actual values of Nn and cos 9, and only weakly dependent on the value of \. 

One considerable uncertainty in our analysis concerns the UV absorption properties of nano-Si material. 
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T (K) 

Fig. 9. — Temperature distribution functions for Si/Si02 spheres of various radii, with the Si cores assumed to be (a) neutral 
and (b) charged with Z = +1, for grains illuminated by the local interstellar radiation field (Mathis, Mezger, & Panagia 1983). 
Curves are labelled by the Si core radius agi 0.8a where a is the Si/SiC>2 radius). For agi =10, 20, 30, 40A a dot indicates 
the first excited state (see Draine & Li 2001). 

In the literature, it is often stated that, as a result of quantum confinement, the optical-UV absorption of 
nano-Si is substantially reduced in comparison with those of bulk-Si. If this is true, the SNP model would 
require an even larger Si abundance to account for the ERE, and the limitations placed on this model would 
be even more severe. 

Another uncertainty concerns the IR emission properties of nano-Si crystals and nano-SiC>2, which have 
not yet been experimentally investigated. In the present work, our results are mainly based on the IR 
emissivities of bulk Si and bulk SiC>2. Nano-materials are expected to be more IR-active than their bulk 
counterparts (due to the symmetry-breaking surfaces) and therefore the detailed emission spectrum will be 
modified. However, the starlight energy absorbed by the grains has to be reradiated in the infrared, so the 
absence of detected spectral features implies upper limits on the SNP abundance. 

Our analysis of NGC 2023 required the value of the foreground extinction A\ to "deredden" the ERE 
measurements. Unfortunately, A\ has not been directly measured at the location where the ERE observations 
have been made. It would be extremely valuable to have ERE measurements made at a position where a 
reliable extinction has been determined - such as the emission bar 80" S of HD 39037 - or to have both IR 
and far-red observations of H2 emission (allowing A\ to be determined) at a location where the ERE has 
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Fig. 10. — As in Fig. 8, but for Si/Si0 2 SNPs with Z Si = 15 ppm. 

been measured. 

Our modelling has assumed isolated "free-flying" SNPs; we have demonstrated that existing DIRBE 
photometry for the diffuse ISM rules out the hypothesis that the observed ERE is due to such particles. Since 
isolated SNPs are ruled out, it is important to note that in the laboratory, photoluminescence is observed for 
SNPs located on a substrate (Ledoux et al. 2000). An SNP located on the surface of a large interstellar grain 
would have its temperature - and therefore its infrared emission - determined by the temperature of the 
"host" grain. In the case of NGC 2023, the radiation field is so intense that the "host" temperature would 
probably be similar to the temperature of the Si/SiC>2 SNPs, so the 20/im SiC>2 emission feature should still 
be prominent, though its strength could be increased or decreased, depending on the optical properties of 
the "host" grain. If the SNPs are pure Si, however, the 16.4/mi emissivity would be reduced if the SNPs 
are attached to a host, since the host grain would not be as hot as the Si SNPs, which are poor thermal 
emitters. In the diffuse ISM, the 16.4 and 20/xm emission would be completely suppressed if the SNPs are 
attached to larger host grains, with temperatures Tra 15 — 20 K. Thus SNPs attached to larger host grains 
may be a viable explanation for the observed ERE. 

The photoexcitation rate 7 for SNPs on the surfaces of a > 100 A grains would probably be lower 
than for free-flying SNPs, which would in turn increase the abundance of Si in SNPs required to explain 
the observed ERE. Since 42% or more of [Si/H] Q is required even for free-flying SNPs, the Si abundance 
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Fig. 11. — Grain shape effects on the model IR emission spectra calculated from spherical grains (thin solid line), 5:1 prolate 
grains (dotted line), 5:1 oblate grains (short-dashed line), and grains with a distribution of spheroidal shapes (long-dashed line) 
illuminated by the NGC 2023 radiation field. All grains are taken to (1) be neutral; (2) have an equivalent sphere-volume Si 
core size a~* = 17.5A ( we choose 17. 5A because Smith & Witt [2001] found that the observed ERE in the diffuse ISM can be 
best explained by SNPs with asi = 17. 5A); (3) have the Si core taking over 51% of the total grain volume. Also plotted are the 
ISO SWS spectrum (heavy solid line; Verstraete et al. 2001) and the KAO photometry (stars; Harvey et al. 1980). It is evident 
that the grain shape effect does not affect our conclusions. 

demands would be difficult to accommodate if the SNPs are attached to a > 100 A grains. 



5. Conclusions 

Table 2 summarizes our results for the abundances of SNPs: upper limits based on nondetection of 16.4 
and 20/im features, and lower limits (since fypL < 1) if neutral SNPs are to account for the observed ERE. 
Note charged SNPs are not able to emit ERE (Smith & Witt 2001). 

For pure (both neutral and charged) Si SNPs, nondetection of a 16.4/zm emission feature in NGC 2023 
implies ?7pl > 0.017 x l0 a4A ° ^». Thus for A).68 M m = 1-2, pure neutral Si SNPs could account for the 
observed ERE in NGC 2023 provided t] PL > 5%, but if Ao.esftm = 3.2, i] PL > 32% would be required to 
account for the ERE. We note in the laboratory high photoluminescence efficiencies are seen only when the 
Si surfaces are "passivated" by an oxide coating. 

An oxide coating, however, will introduce a strong emission feature at 20/im which is not seen in 
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Table 1: ERE Emissivitics of SNPs 
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neutral Si 
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7.43 x 10- 6 (x/2000) 


4.94 x 10- 9 (x/l-23) 




imod ( Cr S S_1 Sr ~ 


-1 H - 


-iy 


1.70 x 10- 18 77 PL Zsi 


1.32 x 10- 21 77 PL Z Si 




(Mabs (eV) 






9.30 


9.29 


charged Si 


7Si (s- 1 Si" 1 ) 






7.04 x 10- 6 (x/2000) 


4.68 x 10- 9 (x/l-23) 






1 H" 


a ) 


1.62 x 10- 18 77 PL Zsi 


1.25 x 10- 21 77 PL Zsi 




(hv) a b s (cV) 






8.59 


8.59 


neutral Si/SiC>2 


7Si (s- 1 sr 1 ) 






6.19 x 10- 6 (x/2000) 


3.87 x 10- 9 (x/l-23) 




Jmod ( Cr S S_1 Sr ~ 


-1 H - 




1.42 x 10- 18 77 PL Zsi 


1.04 x 10- 21 77 PL Zsi 




(Mabs (cV) 






8.65 


8.62 


charged-Si/Si0 2 


7si (s- 1 Si" 1 ) 






6.00 x 10- 6 (x/2000) 


3.77 x 10- 9 (x/l-23) 




imod ( Cr § S_1 Sr ~ 


-1 H - 


- 1 ) 


1.37 x 10- 18 77 PL Zsi 


1.00 x 10- 21 77 PL Z Si 



"For NGC 2023, j™ E is estimated assuming a UV-illuminated slab with Nu = 2 X 10 21 cm 2 with limb-brightening 1/ cos 6 = 2 
(see text). 

b (hv) a b s = J Jj0 ? >1 Cah s {hv)cu v dv I J 550 » A C ^ s {hv)(cu v I hv)dv where cu„ is the starlight intensity, 
yizri y 

c 7si = ^sl 1 C ab .(fti/)(cu„/fti/)di/. 

d Jmod = < ftly >ERE7Si»yPLZsi/4vr. 



Table 2: IR Limits on Si/H (ppm) in Silicon Nanoparticles 







neutral Si 


charged Si a 


Region 


Composition 


Si/H (ppm) from Si/H (ppm) from 
observed IR emission observed ERE intensity 15 


Si/H (ppm) from Si/H (ppm) from 
observed IR emission observed ERE intensity 13 


NGC 2023 


Pure Si 


< 0.8 0.014 r)~* x 10°' 4A ° 


< 1.5 0.014 r t pl x 10°' 4A o.68 M m 

< 0.2 0.017 r)pl x 10 0<L4 o.68,um 




Si/Si0 2 


< 0.2 0.016 »7pL X 10°' 4A 0.68 M m 


Diffuse ISM 


Pure Si 


< 4 12 ?jpV 


< 12 12 TjpY. 




Si/SiQ 2 


< 2 15 r^ 1 


< 2 15 



"Note charged SNPs do not luminesce (i.e. tjpl = 0). 

Vl < 1 is the photoluminescence efficiency. Lower limits are obtained by setting ijpl 
foreground extinction at 0.68/im (see text). 



= 1- A . 68Mm is the 60"S NGC 2023 
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NGC 2023. For Si/SiC>2 nanoparticles to explain the ERE, nondetection of the 20/im feature requires a 
luminescence efficiency fypL > 0.081 x 10°- 4j4 ° f^™. Photoluminescence efficiencies as high as 50% have been 
reported for Si/SiC>2 nanoparticles (Wilson, Szajowski, & Brus 1993), and in principle ?7pl could approach 
100%. However, if A .e8^m ^ 2.8, the ERE could not be due to Si/Si02 nanoparticles since the required tjpl 
would then exceed 100%. 

The ERE emissivity of the diffuse ISM (Gordon et al. 1998; Szomoru & Guhathakurta 1998) requires 
> 42% [Si/H] Q in neutral Si/Si02 SNPs, even for photoluminescence efficiency ?]pl — > 100%. Such a high 
abundance of SNPs in the diffuse ISM is difficult to reconcile with the evidence that a substantial fraction of 
interstellar Si is in fact locked up in amorphous silicate grains (e.g., see Weingartner & Draine 2001). This 
difficulty is exacerbated if interstellar abundances are significantly subsolar, as has been argued (see Snow 
& Witt 1996). 

We have calculated the IR emission which such particles would produce if they are free-flying. Existing 
DIRBE 25^tm photometry appears to already rule out such high abundances of SNPs, indicating that if SNPs 
are responsible for the ERE from the diffuse ISM, they must either be in a > 50 A clusters, or attached 
to larger "host" grains. Future observations by SIRTF will be even more sensitive to the abundance of 
free-flying SNPs in the diffuse interstellar medium. 

In light of the controversy over the optical properties of nano-Si material, we call for further experimental 
studies to determine the dielectric functions of nano-Si material either directly from nano-Si crystals, or from 
porous-Si with the effects of voids and surface oxidation well determined and removed. 

Finally, we also note that, if the surface dangling bonds of Si grains are passivated by H atoms and 
if they are present in interstellar space in quantities sufficient to account for the observed ERE, we would 
expect to see the Si-H 15.6^m wagging, 11.6/xm bending, and even 5/xm stretching bands (Adachi 1999). 

We especially thank A.N. Witt for invaluable comments and helpful suggestions, and L. Verstraete for 
providing us with the ISO spectrum of NGC 2023 in advance of publication. We also thank S. Adachi, 
N. Koshida, P. Maddalena, and V.G. Zubko for sending us their optical constants of bulk Si and/or nano- 
silicon; S. Adachi, L.J. Allamandola, J.M. Greenberg, Th. Hcnning, N. Koshida, D. Kovalev, P. Maddalena, 
L. Pavesi, and K. Sellgren for helpful suggestions; and R.H. Lupton for the availability of the SM plotting 
package. This research was supported in part by NASA grant NAG5-7030 and NSF grants AST-9619429 
and AST-9988126. 
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